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Abstract. 

The attempts to find the right (underlying) theory for the nuclear force have 
a long and stimulating history. Already in 1953, Hans Bethe stated that 
"more man-hours have been given to this problem than to any other scien- 
tific question in the history of mankind". In search for the nature of the 
nuclear force, the idea of sub-nuclear particles was created which, eventu- 
ally, generated the field of particle physics. I will review this productive 
history of hope, error, and desperation. Finally, I will discuss recent ideas 
which apply the concept of an effective field theory to low-energy QCD. 
There are indications that this concept may provide the right framework to 
properly understand nuclear forces. 



1 Historical perspective 

The theory of nuclear forces has a long history (cf. Table 1). Based upon the 
seminal idea by Yukawa [1], first field-theoretic attempts to derive the nucleon- 
nucleon (NN) interaction focused on pion-exchange. While the one-pion ex- 
change turned out to be very useful in explaining NN scattering data and the 
properties of the deuteron [2], multi-pion exchange was beset with serious am- 
biguities [3,4]. Thus, the "pion theories" of the 1950s are generally judged 
as failures — for reasons we understand today: pion dynamics is constrained by 
chiral symmetry, a crucial point that was unknown in the 1950s. 

Historically, the experimental discovery of heavy mesons [5] in the early 
1960s saved the situation. The one-boson-exchange (OBE) model [6,7] emerged 
which is still the most economical and quantitative phenomenology for describ- 
ing the nuclear force [8,9]. The weak point of this model, however, is the scalar- 
isoscalar "sigma" or "epsilon" boson, for which the empirical evidence remains 
controversial. Since this boson is associated with the coiTelated (or resonant) ex- 
change of two pions, a vast theoretical effort that occupied more than a decade 
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Table 1. Seven Decades of Struggle: The Theory of Nuclear Forces 
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and beyond 
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was launched to derive the 27r-exchange contribution to the nuclear force, which 
creates the intermediate range attraction. For this, dispersion theory as well as 
field theory were invoked producing the Paris [10, 11] and the Bonn [7, 12] po- 
tentials. 

The nuclear force problem appeared to be solved; however, with the discov- 
ery of quantum chromo-dynamics (QCD), all "meson theories" were relegated 
to models and the attempts to derive the nuclear force started all over again. 

The problem with a derivation from QCD is that this theory is non- 
perturbative in the low-energy regime characteristic of nuclear physics, which 
makes direct solutions impossible. Therefore, during the first round of new at- 
tempts, QCD-inspired quark models [13| became popular. These models are 
able to reproduce qualitatively and, in some cases, semi-quantitatively the gross 
features of the nuclear force [14, 15]. However, on a critical note, it has been 
pointed out that these quark-based approaches are nothing but another set of 
models and, thus, do not represent any fundamental progress. Equally well, one 
may then stay with the simpler and much more quantitative meson models. 

A major breakthrough occurred when the concept of an effective field theory 
(EFT) was introduced and applied to low-energy QCD [16]. 

Note that the QCD Lagrangian for massless up and down quarks is chirally 
symmetric, i. e., it is invariant under global flavor SU{2) l x SU{2)ii equivalent 
to SU {2)v X SU {2) A (vector and axial vector) transformations. The axial sym- 
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metry is spontaneously broken as evidenced in the absence of parity doublets 
in the low-mass hadron spectrum. This implies the existence of three massless 
Goldstone bosons which are identified with the three pions (tt^, tt"). The non- 
zero, but small, pion mass is a consequence of the fact that the up and down 
quark masses are not exactly zero either (some small, but explicit symmetry 
breaking). Thus, we arrive at a low-energy scenario that consists of pions and 
nucleons interacting via a force governed by spontaneously broken approximate 
chiral symmetry. 

To create an effective field theory describing this scenario, one has to write 
down the most general Lagrangian consistent with the assumed symmetry prin- 
ciples, particularly the (broken) chiral symmetry of QCD [16]. At low energy, 
the effective degrees of freedom are pions and nucleons rather than quarks and 
gluons; heavy mesons and nucleon resonances are "integrated out". So, the cir- 
cle of history is closing and we are back to Yukawa's meson theory, except that 
we have learned to add one important refinement to the theory: broken chiral 
symmetry is a crucial constraint that generates and controls the dynamics and 
establishes a clear connection with the underlying theory, QCD. 

It is the purpose of the remainder of this paper to describe the EFT approach 
to nuclear forces in more detail. 

2 Chiral perturbation theory and the hierarchy of nuclear forces 

The chiral effective Lagrangian is given by an infinite series of terms with in- 
creasing number of derivatives and/or nucleon fields, with the dependence of 
each term on the pion field prescribed by the rules of broken chiral symmetry. 
Applying this Lagrangian to NN scattering generates an unlimited number of 
Feynman diagrams. However, Weinberg showed [17] that a systematic expan- 
sion exists in terms of {Q/A^Y, where Q denotes a momentum or pion mass, 
A-^ K 1 GeV is the chiral symmetry breaking scale, and u > (cf. Figure 1). 
This has become known as chiral perturbation theory (xPT). For a given order 
the number of terms is finite and calculable; these terms are uniquely defined and 
the prediction at each order is model-independent. By going to higher orders, 
the amplitude can be calculated to any desired accuracy. 

Following the first initiative by Weinberg [17], pioneering work was per- 
formed by Ordonez, Ray, and van Kolck [18, 19] who constructed a NN potential 
in coordinate space based upon xPT at next-to-next-to-leading order (NNLO; 
f = 3). The results were encouraging and many researchers became attracted to 
the new field. Kaiser, Brockmann, and Weise [20] presented the first model- 
independent prediction for the NN amplitudes of peripheral partial waves at 
NNLO. Epelbaum et al. [21] developed the first momentum-space NN potential 
at NNLO, and Entem and Machleidt [22] presented the first potential at N^LO 
{u = 4). 

In xPT, the NN amplitude is uniquely determined by two classes of contri- 
butions: contact terms and pion-exchange diagrams. There are two contacts of 
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2N Force 3N Force 4N Force 




Figure 1. Hierarchy of nuclear forces in xPT- Solid lines represent nucleons and dashed 
lines pions. Further explanations are given in the text. 

order [0{Q'^)] represented by the four-nucleon graph with a small-dot ver- 
tex shown in the first row of Figure 1 . The corresponding graph in the second 
row, four nucleon legs and a solid square, represent the seven contact terms of 
0{Q^). Finally, at 0{Q^), we have 15 contact contributions represented by a 
four-nucleon graph with a solid diamond. 

Now, turning to the pion contributions: At leading order [LO, 0{Q^), v = 
0], there is only the welUcnown static one-pion exchange, second diagram in the 
first row of Figure 1 . Two-pion exchange (TPE) starts at next-to-leading order 
(NLO, I' = 2) and all diagrams of this leading-order two-pion exchange are 
shown. Further TPE contributions occur in any higher order. Of this sub-leading 
TPE, we show only two representative diagrams at NNLO and three diagrams 
at N^LO. The TPE at N^LO has been calculated first by Kaiser [23]. All 27r 
exchange diagrams/contributions up to N'^LO are summarized in a pedagogical 
and systematic fashion in Ref. [24] where the model-independent results for NN 
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scattering in peripheral partial waves are also shown. 

Finally, there is also three-pion exchange, which shows up for the first time 
at N^LO (two loops; one representative Stt diagram is included in Figure 1). In 
Ref. [25], it was demonstrated that the 3tt contribution at this order is negligible. 

One important advantage of xPT is that it makes specific predictions also for 
many-body forces. For a given order of xPT, two-nucleon forces (2NF), three- 
nucleon forces (3NF), ... are generated on the same footing (cf. Figure 1). At 
LO, there are no 3NF, and at next-to-leading order (NLO), all 3NF terms can- 
cel [17, 26]. However, at NNLO and higher orders, well-defined, nonvanishing 
3NF occur [26, 27]. Since 3NF show up for the first time at NNLO, they are 
weak. Four-nucleon forces (4NF) occur first at N'^LO and, therefore, they are 
even weaker. 



3 Chiral NN potentials 



The two-nucleon system is non-perturbative as evidenced by the presence of 
shallow bound states and large scattering lengths. Weinberg [17] showed that 
the strong enhancement of the scattering amplitude arises from purely nucle- 
onic intermediate states. He therefore suggested to use perturbation theory to 
calculate the NN potential and to apply this potential in a scattering equation 
(Lippmann-Sch winger or Schrodinger equation) to obtain the NN amplitude. 
We follow this philosophy. 

Chiral perturbation theory is a low-momentum expansion. It is valid only 
for momenta Q <C w 1 GeV. Therefore, when a potential is constructed, 
all expressions (contacts and irreducible pion exchanges) are multipUed with a 
regulator function. 



exp 



2n 



2n 



(1) 



where p and p' denote, respectively, the magnitudes of the initial and final nu- 
cleon momenta in the center-of-mass frame; and A <c A^. The exponent 2n 



Table 2. x^^datum for the reproduction of the 1999 np database below 290 MeV by 
various np potentials. (A = 500 MeV in all chiral potentials.) 



Bin (MeV) 


# of data 


N-'LO" 


NNLO" 


NLO" 


AVI 8= 


0-100 


1058 


1.06 


1.71 


5.20 


0.95 


100-190 


501 


1.08 


12.9 


49.3 


1.10 


190-290 


843 


1.15 


19.2 


68.3 


1.11 


0-290 


2402 


1.10 


10.1 


36.2 


1.04 



"Reference [22]. 
''Reference [28]. 
"^Reference [29]. 
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is to be chosen such that the regulator generates powers which are beyond the 
order at which the calculation is conducted. 

NN potentials based upon xPT at NNLO [21,28] are poor in quantitative 
terms; they reproduce the NN data below 290 MeV lab. energy with a x^/datum 
of more than 20 (cf. Tables 2 and 3). As shown first by Entem and Machleidt 
in 2003 [22], one has to go to order N^LO to obtain a NN potential of accept- 
able accuracy. For a more recent construction of an N^LO NN potential, see 
Ref. [31]. 

For an accurate fit of the low-energy pp and np data, charge-dependence 
is important. Charge-dependence up to next-to-leading order of the isospin- 
violation scheme (NL0, in the notation of Ref. [32]) includes: the pion mass 
difference in OPE and the Coulomb potential in pp scattering, which takes care 
of the L0 contributions. At order NL0 we have pion mass difference in the 
NLO part of TPE, exchange [33], and two charge-dependent contact interac- 
tions of order which make possible an accurate fit of the three different 
scattering lengths, app, Unn, and a„p. 

In the optimization procedure, we fit first phase shifts, and then we refine 
the fit by mininoizing the obtained from a direct comparison with the data. 
The /datum for the fit of the np data below 290 MeV is shown in Table 2, 
and the corresponding one for pp is given in Table 3. The x^ tables show the 
quantitative improvement of the NN interaction order by order in a dramatic 
way. Even though there is considerable improvement when going from NLO to 
NNLO, it is clearly seen that N-^LO is needed to achieve an accuracy comparable 
to the phenomenological high-precision Argonne Vis potential [29]. Note that 
proton-proton data have, in general, smaller errors than np data which expliiins 
why the pp x^ are always larger. 

The phase shifts for np scattering below 300 MeV lab. energy are displayed 
in Figure 2. What the x^ tables revealed, can be seen graphically in this figure. 
The phase shifts are a particularly good example; NLO (dotted line) is 
clearly poor. NNLO (dash-dotted line) brings improvement and describes the 
data up to about 100 MeV. The difference between the NLO and NNLO curves is 
representative for the uncertainty at NLO and, similarly, the difference between 



Table 3. x^/datum for the reproduction of the 1999 pp database below 290 MeV by 
various pp potentials. (A = 500 MeV in all chiral potentials.) 



Bin (MeV) 


# of data 


N^LO" 


NNLO" 


NLO" 


AV18'= 


0-100 


795 


1.05 


6.66 


57.8 


0.96 


100-190 


411 


1.50 


28.3 


62.0 


1.31 


190-290 


851 


1.93 


66.8 


111.6 


1.82 


0-290 


2057 


1.50 


35.4 


80.1 


1.38 



"Reference [22]. 
''See footnote [30]. 
^Reference [29]. 
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Figure 2. np phase parameters below 300 MeV lab. energy for partial waves with J < 2. 
The thick solid (dashed) line is the result by Entem and Machleidt [22] at N^LO using 
A = 500 MeV (A = 600 MeV). The thin dotted and dash-dotted lines are the phase shifts 
at NLO and NNLO, respectively, as obtained by Epelbaum et al. [28] using A = 500 
MeV. The solid dots show the Nijmegen multienergy np phase shift analysis [34], and 
the open circles are the GWUA^PI single-energy np analysis SM99 [35]. 



NNLO and N^LO reflects the uncertainty at NNLO. Obviously, at N^LO (A = 
500 MeV, thick solid line) we have a good description up to 300 MeV. An idea 
of the uncertainty at N'^LO can be obtained by varying the cutoff parameter A. 
The thick dashed line is N'^LO using A = 600 MeV. In most cases, the latter two 
curves are not distinguishable on the scale of the figures. Noticeable differences 
occur only in ^£>2, ^-Pj, and above 200 MeV. 

4 Chiral three-nucleon forces 

As noted before, an important advantage of the EFT approach is that it creates 
two- and many-body forces on an equal footing. The first non-vanishing 3NF 
terms occur at NNLO and are shown in Figure 1 (row 'Q'^/NNLO', column 
'3N Force'). There are three diagrams: the TPE, OPE, and 3N-contact interac- 
tions [27]. The TPE 3N-potential is given by 

T/3NF _ f i (^» • qi}{^j ■ Qj) pa/3 „8 
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with (fi =pi~pi, where pi and pi ' are the initial and final momenta of nucleon i , 
respectively, and 



F«/3^ca/3 _4cim|. 2C3 
ijk " f2 f2 

J TV J TT 



+ ^^e«'57^7^,.[^.X5,]. (3) 



The vertex involved in this 3NF term is the two-derivative tt-kNN vertex (large 
solid dot in Figure 1) which we encountered already in the TPE contribution to 
the 2N potential at NNLO. Thus, there are no new parameters and the contribu- 
tion is fixed by the LECs used in NN. The OPE contribution is 

and, finally, the 3N contact term reads 

V^r = Elj:^rr.. (5) 

The last two 3NF terms involve two new vertices (that do not occur in the 2N 
problem), namely, the ttNNNN vertex with parameter D and a 6iV vertex with 
parameters E. One way to pin down the two new parameters is to fit them to 
the triton and the ^He binding energies. Once D and E are fixed, the results 
for other 3N, 4N, . . . observables are predictions. Results for 3N scattering ob- 
servables are reported in Refs. [36, 37]. Spectra of Ught nuclei are calculated 
in Refs. [38, 39]. Concerning the famous 'Ay puzzle', the above 3NF terms 
yield some improvement of the predicted nucleon-deuteron analyzing powers, 
however, the problem is not resolved. 

One should note that there are additional 3NF terms at NNLO due to rela- 
tivistic corrections (l/M^ corrections) that have not yet been included in any 
calculation. However, there are all reasons to believe that these contributions 
will be very small, probably negligible. It is more Ukely that the problem with 
the chiral 3NF is analogous to the one with the chiral 2NF: namely, NNLO is 
insufficient and for sufficient accuracy one has to proceed to N^LO. Two 3NF 
topologies at N^LO are indicated in Figure 1. The N^LO 3NF, which does not 
depend on any new parameters, is presently under development. 



5 Conclusions 



The EFT approach to nuclear forces is a modern refinement of Yukawa's meson 
theory. It represents a scheme that has an intimate relationship with QCD and 
allows to calculate nuclear forces to any desired accuracy. Moreover, nuclear 
two- and many-body forces are generated on the same footing. 
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At N^LO [22], the accuracy is achieved that is necessary and sufficient for 
microscopic nuclear structure. First calculations applying the N'^LO NN poten- 
tial in the no-core shell model [40^2], the coupled cluster formaUsm [43^8], 
and the unitary-model-operator approach [49] have produced promising results. 

The 3NF at NNLO is known [27] and has had first successful applications 
in few-nucleon reactions [36,37] as well as the structure of light nuclei [38,39]. 
The 3NF at N^LO is under construction. 

It may be too early to claim that the nerver-ending story is coming to an end, 
but the story is certainly converging. 
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